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Abstract 

(N 

A linked cluster expansion suitable for the treatment of ground-state properties of complex nuclei, 
as well as of various particle-nucleus scattering processes, has been used to calculate the ground- 
state energy, density and momentum distribution of 16 O and 40 Ca using realistic interactions. 
', First of all, a benchmark calculation for the ground-state energy has been performed using the 

truncated V8' potential, and consisting in the comparison of our results with the ones obtained 
by the Fermi Hypernetted Chain approach, adopting in both cases the same mean field wave 
functions and the same correlation functions. The results exhibited a nice agreement between 
the two methods. Therefore, the approach has been applied to the calculation of the ground- 
' state energy, density and momentum distributions of ie O and 40 Ca using the full V8 1 potential, 

finding again a satisfactory agreement with the results based on more advanced approaches where 
higher order cluster contributions are taken into account. It appears therefore that the cluster 
expansion approach can provide accurate approximations for various diagonal and non diagonal 
density matrices, so that it could be used for a reliable evaluation of nuclear effects in various 
medium and high energy scattering processes off nuclear targets. The developed approach can 
be readily generalized to the treatment of Glauber type final state interaction effects in inclusive, 
semi-inclusive and exclusive processes off nuclei at medium and high energies. 
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I. INTRODUCTION 



The knowledge of the nuclear wave function, in particular its most interesting and un- 
known part, viz the correlated one, which is predicted by realistic many-body calculations 
to strongly deviate from a mean field description, is not only a prerequisite for understand- 
ing the details of bound hadronic systems, but is becoming at present a necessary condition 
for a correct description of medium and high energy scattering processes off nuclear targets; 
these, in fact, represent nowadays an efficient tool for the investigation of several high energy 
problems, e.g. color transparency, hadronization, the properties of dense hadronic matter, 
etc. which manifest themselves only in the nuclear medium. The necessity of an accurate 
treatment of the effects of the medium in high energy scattering process is becoming a rel- 
evant issue in hadronic physics. The problem is not trivial, for one has first to solve the 
many body problem and then has to find a way to describe scattering processes in terms 
of realistic many-body wave functions. The difficulty mainly arises because even if a reli- 
able and manageable many-body description of the ground state is developed, the problem 
remains of the calculation of the final state. In the case of few-body systems, a consistent 
treatment of initial state correlations (ISC) and final state interaction (FSI) is nowadays 
possible at low energies by solving the Schrodinger equation for the bound and continuum 
states (see e.g.Q, 0, [| and References therein quoted), but at high energies, when the num- 
ber of partial waves sharply increases and nucleon excitations can occur, the Schrodinger 
approach becomes impractical and other methods have to be employed. Moreover, in the 
case of complex nuclei, additional difficulties arise due to the approximations which are still 
necessary to solve the many-body problem. As a matter of fact, in spite of the relevant 
progress made in recent years in the calculation of the properties of light nuclei (see e.g. 
"\ El y, Q, )j much remains to be done, also in view that the results of very sophisticated 
calculations {e.g. the variational Monte Carlo ones |s|), show that the wave function which 
minimizes the expectation value of the Hamiltonian provides a very poor nuclear density; 
moreover, the structure of the best trial wave function is so complicated, that its use in the 
calculation of various processes at intermediate and high energies appears to be not easy 
task. It is for this reason that the evaluation of nuclear effects in medium and high energy 
scattering processes is usually carried out within simplified models of nuclear structure. As 
a matter of fact, ISC are often introduced by a procedure which has little to recommend 
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itself, namely the expectation value of the transition operator is evaluated with shell model 
(SM) uncorrelated wave functions and the initial two-body wave function describing the 
independent relative motion of two nucleons is replaced by a phenomenological correlated 
wave function. Recently, however, important progress has been made, in that ISC have 
been introduced from the beginning using correlated wave functions and cluster expansion 
techniques. Central Jastrow-type correlations have been often used to investigate the effects 
of ISC on various scattering processes off complex nuclei induced by medium energy leptons 
like, e.g. A(e,e')X 10|, A(e,e'p)X Q and A(e,e'2N)X Q processes. Calculations of in- 
clusive electron scattering has also been performed using realistic many body wave functions 
and spectral functions within various approximations [3,0], an d non central correlations 
have been recently introduced in the calculation of A(e,e'p)B and A(e,e'2N)B in 
spite of this progress, further work remains to be done to achieve a full consistent treat- 
ment of both ISC and FSI in intermediate and high energy scattering off complex nuclei. 
This would be particularly urgent, as far as various high energy phenomena are concerned, 
e.g. exclusive processes at high momentum transfer [171 ] . inclusive |1S| and semi inclusive 
[19j ] hadron production in Deep Inelastic Scattering, and others, which might also require 
a careful treatment of nuclear effects. As a matter of fact, a recent calculation [2^ of the 
integrated nuclear transparency in the processes 16 0(e,e'p)X and 40 Ca(e,e'p)X performed 
within a cluster expansion approach, including realistic central and tensor correlations, shows 
that the results do depend both on the SM and ISC parameters, which have therefore to be 
fixed from firm criteria, e.g. from the calculation of the static properties of nuclei, like the 
ground-state energy and the density distribution. 

For such a reason, we have undertaken the calculation of the basic ground-state proper- 
ties (energies, densities and momentum distributions) of complex nuclei within a framework 
which can easily be generalized to the treatment of various scattering processes, keeping the 
basic features of ISC as predicted by the structure of realistic Nucleon-Nucleon (NN) inter- 
actions. Our approach is presented in detail in this paper, which is organized as follows: in 
Section II some basic ideas concerning the application of the cluster expansion techniques to 
the approximate solution of the nuclear many-body problem are recalled; the cluster expan- 
sion used in the calculations is described in Section III; the ground-state energy calculations 
for 16 and 40 Ca are presented in Section IV, where the results a benchmark calculation 
aimed at a comparison of our results with the results obtained within the Fermion Hypernet- 
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ted Chain (FHNC) approach is presented; the results of the calculations of the charge 

densities and momentum distributions are given in Sections V and VI, respectively; the dia- 
grammatic representation of the latter quantities within the cluster-expansion approach are 
illustrated in Section VII; the Summary and Conclusions are presented in Section VIII. 
Preliminary results of our calculations have been presented in Ref . j^li I22I ] . 

II. THE CORRELATED WAVE FUNCTIONS 

It is well known that if nuclei are considered to be aggregates of point-like nucleons with 
the same properties and interactions as the free ones, and, moreover, all degrees of freedom 
but nucleonic ones are frozen, the nuclear many-body problem reduces to the search of the 
eigenvalues and eigenfunctions of the following Schrodinger equation 



A 2 

' Vn(l,2,...,A) = E n i/; n (l,2,...,A) (1) 



i=i 



2M N 



where M^r is the nucleon mass and the effective interaction includes many-body interactions 
between the costituents, i.e. 

V eff (l,2,...,A) = 5> 2 (*,j) + E vs(i,j,k) + ... + v A (l,2,...,A). (2) 

i<j i<j' 



Within the so called standard model of nuclei [24(, which will be considered from now on, 
many-body interactions are disregarded and Eq. is solved keeping only the two-body 
interaction $2(1, j), whose form is determined from two-body bound and scattering data. 
We will consider therefore the following nuclear Hamiltonian: 

H = f + V = E V 2 + E , (3) 

ZM N i= i i<j 

where the vector x denotes the set of nucleonic degrees of freedom, i.e. x = (r; cr; r), with 
r, a and r denoting the spatial, spin and isospin coordinates, respectively. We will try to 
find the solution of the Schrodinger equation pertaining to the ground-state of the nucleus, 
i.e.: 

HiPo = E tfj (4) 

and to this end we will look for the ground-state wave function (WF) tp Q which minimizes 
the expectation value of the Hamiltonian 

<*> - ™ * K. (5) 



As trial WF we will use a correlated WF of the following form 

i/) (xx,...,xa) = F(x l ,...,x A )<p (x 1 ,...,x A ) , (6) 

where <p Q is a SM , mean-field WF describing the independent particle motion, and F is a 
symmetrized correlation operator, which generates correlations into the mean field WF; the 
correct symmetry of the WF is guaranteed by <ft . 

As in any variational approach, the central problem here is to give an explicit form to 
the trial WF (Eq. ©); whereas for <ft any realistic SM WF can be considered a physically 
sound approximation, the choice of the form of the operator F is not clear apriori. However 
one can be guided by the knowledge of the basic features of the force acting between the 
considered hadrons. Nowadays the nucleon-nucleon interaction can be cast in the following 
form ( 3): 

v = x> (n) K-)0|;\ (7) 

n=l 

where = — rA is the relative distance of nucleons % and j, and n, ranging up to N = 18, 



labels the state-dependent operator O 



(n). 



(n) 



1 , <Tj • a-j , Sij , (S ■ L)^ , L 2 , L 2 cr j • a-j , (S ■ L)\. 



® [i , n ■ tj] (8) 



Accordingly, the operator F is written as 



A 

F(xi, x 2 ... x A ) = Sn/>ii) (9) 

i<j 

with 

N 

/(%) = E / (B) (r y ) /W(r«) = / (n) (^) Ojf • (10) 

n=l 

The variational principle requires the full evaluation of Eq. (j3J) which, obviously, is no easy 
task due to the structure of ipo- We will evaluate the expectation value of the Hamiltonian 



using the cluster expansion techniques 
described in the next sections. 



adopting a specific cluster expansion to be 



III. THE CLUSTER EXPANSION 



The evaluation of the expectation value of H is object of intensive activity which in the 
last few years has produced considerable results: the approximate solution of the Schrodinger 
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equation by means of Monte Carlo methods, for example, has reached a great level of 
accuracy, and the ground-state properties of nuclei with A = 16 have been obtained with a 



full evaluation of Eq. @ [4, 
FHNC approximation p, 



n 



H; exhaustive calculations have also been performed within the 
|9|. Nevertheless, the level of complexity of these calculations 
often prevents the WF to be used with reasonable ease in other nuclear-related problems, 
such as nuclear reactions. Our goal is to present a more economical, but effective method 
for the calculation of the expectation value of any quantum mechanical operator O in the 
many-body ground-state described by the WF ip 01 i.e.: 



(O) 



(ii) 



(V'olV'o) 

with ip having the structure of Eq. © . In the present section we are going to introduce a 
cluster expansion technique in order to evaluate Eq. (fTTj) . To begin with, a generic operator 
O in Eq. (fTTj) will be considered in the following Section, while in the next Section, it will 
be specialized to the Hamiltonian and to the one- and two-body density operators. 

Various types of cluster expansions have been used in the past to calculate the ground- 
state properties of nuclei (see e.g. j^, in these calculations , mainly aimed at inves- 
tigating the convergence of the expansion, simple models of the NN interaction have been 
usually used . In this paper we use an expansion which has never been used previously to 
calculate gropund- state properties of nuclei in terms of realistic interactions. The expansion 



we are going to use has been originally developed in Ref. Q (see also Q and |3jj); the 
main feature of such an expansion is that it is linked and number conserving. The latter 
property means that the normalization of any observable is provided by the normalization 
of the mean field WF, i.e., by the first term of the expansion: the contribution of all other 
terms to the normalization vanishes analytically order by order. The expansion, to be called 
the rj -expansion, has been originally used to obtain the lowest order contribution to the 

i) 13; an d to the one-body mixed density 
30], using central correlations only (i.e., Eqs. © and (fTUj) 
, the lowest order expansion of the OBMD operator has 



diagonal one-body density (OBD) matrix, p^>(r 
(OBMD) matrix, pW^ri) 
with N — 1). Subsequently 



32 

been generalized to take into account also the non central spin-isospin and tensor-isospin 
correlations f 4 and f 6 in Eq. ([TOjl . which turned out to be the most relevant non central 
correlation functions in Nuclear Matter ([33]), as well as all correlations up to N = 6 Q] 1 . 



The approximation which includes only the components n = {1,4,6} is usually referred to as the f$ 
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During the last few years, lowest order expansions have also been applied, within the central 

n 

correlation approximation, to the calculation of the two-body density matrix [3J| and of var- 
ious transition matrix elements appearing in inclusive, A(e, e')X, and exclusive, A(e, e'p)B 
and A(e, e'2N)B, processes (see e.g. Refs. 11> 12 1). As already mentioned, the expan- 



sion has also been used, within the f% approximation, to calculate the nuclear transparency 



in the semi-inclusive process A(e,e'p)X |20j |. To our knowledge, the ^-expansion has never 
been used to calculate the ground-state energy of complex nuclei with a realistic interac- 
tion. It is precisely the central aim of our work to present a detailed report of the results of 
the calculation of the ground-state energy, density and momentum distributions of complex 
nuclei using the ^-expansion and realistic interactions. 

Let us first of all recall the basic features of the expansion. Following the formal expression 
for ip in and taking the correlation operator F as in^J one writes 

f 2 = n/ 2 (%) =n(i+^)) = 

i<j i<j 

= l + J2Vij+ Yl VijVki + ■■■ (12) 

i<j (ij)<(kl) 

where 

m = /?•-! (13) 

and (\r]\ 2 ) will play the role of a small expansion parameter, in that its expectation value 
on the reference state <fi is small. We use the notation rjij = ^(r^). In what follows, when 
dealing with state-dependent operators, we have to bear in mind that they do not commute 
with each other; moreover, since the same operators appear both in the potential and in the 
WF, it is worth defining the following quantities: 

1) (A) = (0 O I A \ 4> ) , where A is an arbitrary quantity; 

2) fjijd = f tJ df t3 - 6- 

3) fjij fjki O ee f i:j f ki 6 fij fki - fij O fij - f k i 6 fki + O; 

and so on, where O is the operator appearing in Eq. (fTTj) . Let us now perform the expansion 
of the expectation value (fTTj). Keeping in mind the described recipes, the quantity F^F is 



approximation, whereas the approximation which includes all correlations up to N = 6 is referred to as 
the fe approximation. 
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expanded both in the numerator and the denominator and all terms containing the same 
number of functions rjij = r](ij) are collected, obtaining 

(6) = O a + O l + 2 + ... + O n + ... + O a , (14) 

At 2-nd order in rj, one has, explicitly, 

= (6) , (16a) 

01 = (Eva 6) - Co <£%>, (166) 

ij ij 

2 = ( E vm 6) - <E mj 6) <E %) + ( 16c ) 

-o (<£M«> -(E%) 2 ) ; ( 16rf ) 

where the term of order n contains f] (/) up to the n-th (^n-th) power. 

Analyzing the structure of Eqs. (j!6dj) one realizes that there are terms which are due to 
the expansion of the denominator 

— — ~ 1 + x + ... , (17) 
1 — x 

e.g. the 2-nd term in Eq. (|16&[) . 

A nice feature of the 77-expansion shows up at this point. Each of the terms in the 
residual formulae presents some linked and unlinked contributions. What we mean by this 
is quite self-explaining when expressed in terms of standard Mayer diagrams according 
to which in each linked term the involved n < A particle coordinates are connected either 
by an / factor, or by Pauli correlations (see Section VII). It is precisely the expansion of 
the denominator which ensures that only linked terms contribute to the overall expectation 
value, all unlinked terms cancelling out amongst themselves. This feature turns out to be 
very convenient from a computational point of view, for it reduces the number of involved 
linked terms and allows one to obtain a very systematic and general procedure. Eventually 
it should be stressed that because of the non-commutative nature of the operators 
involving at least one common particle index, sets of diagrams involving more than two 
particles appears in the expectation value of the Hamiltonian expression within the cluster 
expansion, already at first order. 
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IV. APPLICATION OF THE ^-EXPANSION TO THE NUCLEI lb AND 40 Ca: 
GENERAL FORMULA AND A BENCHMARK CALCULATION WITH TRUN- 
CATED V8' AND [714 INTERACTIONS 



A. General formulae in terms of density distributions 

Given the two-body interaction as in Eq. (J7J), the expectation value of the Hamiltonian 
can be written in the following way j^]: 

E ° = -Y^j dr i [ V ^ (1) ( r i' r, i)]r 1=n + E / dr ^ ^ (n) (^)P(?)(n,r 2 ) , (18) 

where p^(r*i, r[) and p( 2 '(ri, r 2 ) are the OBMD and the TBD matrices, respectively, which 
are defined as the expectation value of the operators 

Pi(n^'i) = E tin ~ ?i) s « ~ *i) II s ( r i - r 'j) (19) 



and 



i.e. 



where 



and 



p 2 (f 1; f 2 ) = E - f <*fo - f 2), (20) 

i<j 



P^iny,) = (MpifarQM) (21) 
pV\r u r 2 ) = (^ \p2(ri,r 2 )\^ ), (22) 

tpo = ip (xi, ...,x A ) (23) 



^ = ^K,...,^) (24) 

and a summation over spin and isospin variables is implicit in Eqs. (121)) and (J22)) . The 
knowledge of the OBMD and TBD matrices allows one to calculate, besides the ground- 
state energy, other relevant quantities like e.g. the density distribution: 

p(r) = pM(r 1 = r' 1 = r), (25) 

the mean square radius of the distribution: 

(r 2 ) = J drr 2 p(r) (26) 
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and, eventually, the nucleon momentum distribution , i.e. the square of the WF in momen- 
tum space which, by definition, reads as follows: 

n(k) = j^y 3 J d ri dr\e^-^ p^r-) . (27) 

The normalization of the OBD, OBMD and TBD matrices are as follows 

Jdrp(r) = A, (28) 

jdr[pV( ri ,r[) = p«(r-i) (29) 

Jdndr.p^ir^r,) = ; (30) 

with the sequential relation 

/dr 2 p< 2 >(n,r 2 ) = ^yVVi) (31) 

linking p^ 2 \ri,r 2 ) to p( l \ri). Accordingly, the normalization of the nucleon momentum 
distribution is 

J dkn(k) = A (32) 

It is useful at this moment to recall the form of p(r±), p^(ri,r[) and p^(ri,r 2 ) 
predicted by the SM (or mean field) approximation. In this case one has ip Q = 
4> = (Aiy 1 / 2 det{(f ai (xj)}, with the single particle (s.p.) orbitals given by <p a (x) = 
ip a (r) x]J 2 ^y 2 ) where a = {a; a; r} = {n, I, m; a; r}. For closed shell nuclei one obtains: 

PsLifi) = 5]|v?a(^i)| 2 = 4p (ri) (33) 

a 

and 

pS#(ri,ri) = E^(xiWxl) = 4p< 1 >(r 1 ,r' 1 ) (34) 
where the sum over a runs over the occupied SM states below the Fermi level, and 

PoiTi) = 5Zko(ri)| 2 (35) 

a 

and 

^(ri.ri) = E^(n)^W). (36) 
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For the TBD matrix one obtains 

Z a/3 

2 

where p (fi) = ( r ii r i)- When Eq. (J21j) is evaluated with the correlated wave functions 



l -A(Ap { ri )p {r 2 ) - tfVi.raJtfVsfcn)), (37) 



(jHJ within the f$ approximation, at first order of the ^-expansion the following expression 
is obtained: 

P (1) (n,r;) = pg^ri) + ^(rx.ri) + pgWi) , (38) 

with 

= J dr 2 [H D (r u ,r V2 )p^\r u r[) Po (r 2 ) ~ H E (r u ,r V2 )p^\ ril r 2 )p^\r 2l r[^ (39) 

dr 2 dr 3 p^\r 1 ,r 2 )[H D (r 23 )p^(r 2 ,r' 1 )p (r 3 ) - H E (r 23 )p^(r 2 , r 3 )£\r 3 , r[)] (40) 

where = |r, — rj\. The subscripts H and 5, whose meaning will be explained in Section 
IVH| stand for hole and spectator, respectively, and 

H D ^\r ih r„) = f {P \n 3 ) f {q \r kl ) cg$(r y , r„) - C^fa, r„) . (41) 

p,q=l 

Here the subscripts D and E stand for direct and exchange, respectively, and the coefficients 
{ r iji r ki)i C E ' q \rij,r k i), whose explicit expressions are given in Appendix C, result 
from the spin and isospin summation, and their explicit dependence upon the coordinates 
originates from the tensor operator. As for the correlated TBD matrix, this can be obtained 
by multiplying the TBD operator in Eq. ({20} by the operators 0\j of Eq. (|8|); the resulting 
TBD matrix, corresponding to the operator reads as follows: 

P$(n.r a ) = ?$(rx,r a ) + p$(n,r 2 ) + P^{r^r 2 ) + p$(n,n) (42) 



-(<l)<l) 4 S) - 4 n) ) P„(ri, r a ) Po (r 2 , r,)) (43) 
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(p) f(<l) r(p»9) _ I v 

),7> D (n),? 1 X 



x p (n, r m ) p (r 2 , r> {2} ) p G (r 3 , r> {3} ) (44) 



(p) jr(g)r»CPi«) _ t v 

(n),P °(n),-P 1 X 



x Po(n, r-p { i } ) p (r 2 , t> {2} ) p G (r 3 , t> {3} ) (45) 



A(A- 



x p (r u r v{1} ) p (r 2 , r v{2} ) p c (r 3 , r v{3} ) p (r 4 , t> {4} ) (46) 



This expression deserves a few explanations; A^)™, B^f v , C^yp and D[X-p are the 
result of the spin-isospin summations and they are in general function of the coordinates; 
the remaining summations over the spatial quantum numbers are then expressed in terms 
of combinations of OBMD matrices (see Appendix C); the subindex V in these factors 
stands for all possible permutations of the states but the unlinked one and the sub-indexes 
V{i} means the corresponding index resulting from the particular permutations; finally, the 
matrices Kr-g) are proper numerical combination of the spin-isospin operators and are 
defined by the following relation 



Aip) Aii) _ K (r) / 4? x 



Sip) Aii) _ ifir) Air) 

r=l 

Note that even if we are dealing with two-body correlations and interactions we end up with 
three- and four-body operators, due to the fact that, e.g. terms like 0(i, 2 )0(i,3), cannot be 
further reduced. Thus, the first order ^-expansion for the energy gets contributions from up 
to four-body clusters. 

iFrom the definition of the nucleon momentum distribution, i.e. Eq. (J27|) . we can obtain 
the expectation value of the kinetic energy operator as follows: 

^ = lkl dkk2n ^- < 48 » 

and Eq. (fTHJ) . finally becomes 

h 2 

I Y 



Eo= Sr I dkk2 < k ^ + E / rfn^ 2U W(r 12 )pg ) ) (r 1 ,r 2 ), (49) 
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with p[^(ri,r 2 ) given by Eq. This is the final expression which has been used to 

calculate the ground-state energy by the following procedure: we have calculated at the same 
order both n(k) and then by placing them in Eq. (J4*§)) and performing the summation 
over n the ground state energy Eq is obtained. Calculations have been performed with a 
given, fixed form for the correlation functions, and considering as variational parameters, 
the parameters of the mean field wave functions. To begin with, in the next Section the 
results of a benchmark calculation aimed at investigating the convergence of the expansion 
will be presented. 

B. A benchmark calculation for ie O ans 40 Ca: comparison between the yy-expansion 
and the Fermion-Hyper-Netted-Chain / Single Operator Chain (FHNC/SOC) ap- 
proach with truncated V8' and J714 interactions 

In order to investigate the convergence of the ^-expansion, we have performed a bench- 
mark calculation consisting in a comparison of Eq. ()49|) with the energy predicted by the 
FHNC/SOC approach. Namely, we have calculated the ground-state properties of 16 and 



and UU 



interactions, respectively 



40 'C a using the first six components of the V8' 
(these model interactions are usually referred to as the truncated V8' and £714 interactions). 
The results we have obtained are compared with the results obtained with FHNC/SOC using 
the same interaction, the same mean field WF's, and the same correlation functions. The six 
correlation functions used in the calculation for 16 , corresponding to the V8' interaction, 
are shown in Fig. Q and the results of the energy calculation are presented in Tables |l] and 
ITT1 It can be seen that the cluster expansion results are very similar to the ones provided by 
the FHNC/SOC method; particularly worth being mentioned is the almost identical value of 
the mean kinetic energy, which means that the nucleon momentum distributions predicted 
by the two methods are very similar. The results of the calculation for 40 C a, correspond- 
ing to the truncated U14 interaction [3jj| and to the mean field and correlation functions 
shown in Fig. |21 are presented in Table IIHI where they are compared with the results of 
the FHNC/SOC approach of Ref. jsj]. Being the mean field wave functions and correlation 
functions the same in the two calculations, any difference between our results and those of 
Ref. Julias to be ascribed, as in the case of 16 0, to the contributions that are left out in the 
cluster expansion. It can be seen that the difference between the two approaches is larger 
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Ill m Ca than in the 16 case, the largest difference arising from the spin-isospin interaction, 
which as a matter of fact is of longer range in m Ca (cf. Fig. |2j). 

To sum up, it seems that the convergence of the ^-expansion for the ground-state energy 
is a satisfactory one. 



V. APPLICATION OF THE //-EXPANSION TO THE NUCLEI 16 AND 40 Ca: 
THE GROUND-STATE ENERGY, RADIUS AND DENSITIES WITH THE FULL 
V8' INTERACTION 

In Ref. p|, using the full V8' interaction which includes the spin-orbit contributions Vj 
and Vs, several ground-state properties of 16 and 40 C a have been calculated within the 
FHNC/SOC approach, namely the ground state energy and the density and momentum 
distributions. For this reason, we have also calculated the ground-state properties of 16 
and 40 C a using the //-expansion and the correlation functions of Ref. 3] which are shown 
in Fig. |3] and 0] The FHNC/SOC calculation of Ref. [7j was performed within the /g 
approximation. We have also used such an approximation but, unlike Ref. [17], we have 
disregarded the v-j and v% components of the V8' interaction. For such a reason a direct 
comparison of the results for the potential energy is not possible, whereas a comparison of 
the average kinetic energy is fully meaningful. The results of the comparison are presented 
in Tables IVl and [VTI for 16 0, and IVII| and IVHTl for m Ca. The most striking feature of the 
correlation functions obtained in Ref. [Tj is the long tail of the tensor-isospin correlation 
function f 6 , which is expected to affect the convergence of the cluster expansion. As a 
matter of fact, it can be seen that whereas the difference in (T) are of the same order as 
in the benchmark calculation, for (V) the situation is not as good as for (T). The reason 
should probably be ascribed to the long tail in / 6 , and the spin-orbit term in the interaction 
which is dropped in the present calculation. As far as the latter is concerned, we have 
estimated the effect of the inclusion of the angular momentum dependent terms, by using 
the nuclear matter results of Ref. [3], and the discrepancy for the nuclear matter case seems 
to be consistent with the discrepancy we found. 

The results we have obtained deserve the following comments 

i) We have compared our results obtained with the truncated V8' and the f G approxima- 
tion but using the full VS' potential which includes the v 7 and v 8 components. Since in 
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both calculations the same mean field WF and correlations functions have been used, 
the differences between the two results have to be ascribed to the terms left out in the 
cluster expansion. Our estimate of the contribution of the V7 and v 8 , based on nuclear 
matter results, shows that this seems indeed to be the case; 

ii) the average kinetic energy obtained in p] agrees with the one obtained by our approach; 
we will indeed show that the momentum distribution, from which the kinetic energy is 
obtained, (see Eq. (|48p). is in very good agreement with the one obtained in Some 
discrepancies are still present as far as the potential energy is concerned, but obtaining 
a full agreement between the lowest order cluster expansion and the FHNC/SOC 
approaches is illusory. 

iii) the overall value of the ground-state energy obtained in this section is reasonably 
closer to the experimental one (~ 8 MeV per nucleon) and it appears that the r\- 
expansion provides a reasonable WF as far as the expectation value of the Hamiltonian 
is concerned. 

By letting V\ = r[ = r in Eq. (|38|) . the matter density at first order of the //-expansion, 
is obtained, i.e. 

p(r) = 4p (r) 

+ J dr 2 [H D {r 12 ) Po {r) Po {r 2 ) - H E {r 12 ) p£\r,r 2 ) p£\r 2 ,r) 
- J dr 2 dr 3 pM(r,r 2 )[H D (r 23 ) Po (r 2 ,r) Po (r 3 ) - H E (r 23 ) p^(r 2 , r 3 ) p?(r 3 , r)} . (50) 



re charge density of the proton 
). Using the mean field WF 



The charge density is obtained by convoluting p(r) with t 
and by correcting for the center-of-mass motion (see e.g. 
and the correlation functions obtained from the ground-state energy calculation with the 
full V8' interaction (cf. Figs. 3 and 4 and Tables V-VIII), the densities shown in Figs. Eland 
El for 16 and for 40 Ca, respectively, are obtained. The results presented in Fig. El and Fig. 
El clearly show that the charge density calculated within the first order //-expansion agrees 
very well with the results obtained in [7] within the FHNC/SOC approach, which indicates 
a very good convergence of the //-expansion as far as the density is concerned. 

It should however be pointed out, that, as first found in [jj, the density calculated with 
mean field WF which minimizes the ground-state energy, strongly disagree with the experi- 
mental density. To cure such a problem, following Ref . (jj , we have recalculated the density 
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varying the mean field parameters to obtain an agreement with the experimental density. 
We take advantage of the fact that, as shown in Fig. [7[ the energy minimum calculated 
within the ^-expansion is a rather shallow one. The results are shown in Figs. [H] and El and 
the comparison with the results of Ref. jx] demonstrate once again the good convergence of 
the ^-expansion. 

The six different two-body densities distributions (Eqs. (42)-(46)) corresponding to the 
first six correlation operators, are shown in Fig. El for 16 (top) and 40 Ca (bottom): each 
of these densities couples with the corresponding component of the realistic potential to 
give the potential energy expectation value. Note that the quantities shown in Fig. ITUlare 
integrated over the center of mass variable , i.e. 

P(?)W =4tt / dRp^ ) (R= l -(r 1 + r 2 ),r=\r\ = \r 1 -r 2 \) . (51) 

VI. APPLICATION OF THE ry-EXPANSION TO THE NUCLEI 16 AND 40 Ca: 
THE NUCLEON MOMENTUM DISTRIBUTIONS 

Using the correlation functions shown in Figs. H3 and EJ and the mean field WF's corre- 
sponding to the best densities shown in Figs. [HI and El , we have calculated the momentum 
distributions given by Eq. (J27j) . with the OBMD matrix p^(ri, r[) given by Eq. (J38|) . The 
results, obtained at first order in rj (the convergence will be discussed later on), are presented 
in Figs. ITTl IT21 [Till and [TH In Figs. ITT1andlT2l our results are compared with the the results 
obtained in Ref. where the same interaction and the same correlation functions have 
)een used, and in case of 16 0, also with the results of Ref. where the AV1A interaction 
Sfij ] and the Variational Monte Carlo approach have been used. These comparisons show 
that: 



. Q and 



1. our results nicely agree with the ones of Refs 

2. short range central correlations do not produce enough high momentum components, 
although they appreciably affect the momentum distributions at k > 2 /m _1 ;the 
inclusion of the tensor operators greatly enhances the high-momentum tail of the 
distribution in the region k > 2 /m _1 . 
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3. the largest effect from non central correlations comes from tensor (Sy) and tensor- 
isospin (Sij Tj ■ Tj) correlations (n = 4 and n = 6 in Eq. (jHj)), the other components 
playing a minor role; thus, as shown in Fig. the f% approximation appears to be a 
rather good one for the calculation of the momentum distributions; 

iv) the satisfactory agreement of our results with the ones of Refs. Q and shows that 
the convergence of the ^-expansion for p(r, r') is a very good one. As a matter of fact, 
we have explicitly evaluated the next order cluster contribution for 16 0; the results, 
reported in Fig. using the f% approximation for the correlations functions, show a 
very good convergence indeed. 

We would like to stress that the good convergence of the momentum distributions is a 
proof of the good convergence of (T) . 



VII. EFFECTS OF CORRELATIONS ON THE CHARGE DENSITY AND MO- 
MENTUM DISTRIBUTIONS: THE DIAGRAMMATIC DESCRIPTION 

Within any body approach based upon the correlated wave function (jUJ), the interaction 
(J7J) and the cluster expansion technique, any quantity can be described by a diagrammatic 
representation, which provides a meaningful definition of correlations and the extent to which 
they affect the given quantity. In particular, the density and the momentum distributions 
are associated to diagrams according to the following rules (see e.g. Ref. 

1. an open dot with the index i denotes the coordinate rf, 

2. a full dot with index i stands for integration over 7*,; 

3. an oriented line joining two dots with indexes i and j denotes p^ir^ rj); 

4. a line beginning from and ending in a dot with index i denotes p (rj); 

5. a dashed line joining two full dots with indexes i and j denotes H D( - E \rij); 

6. two dashed lines joining the open dots with index i and j with the full dot denotes 
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The diagrammatic representations of the diagonal (Eq. (|5Uj) ) and non diagonal (Eq. ((58)1) 
one body density matrices are shown in Figs. EH and El respectively. The meaning of 
hole (H) and spectator (S) contributions introduced in Section II V Al (cf. Eqns. (|39|) and 
()4Ujl )becomes now clear: the first, represented by diagrams 156 (166) describes the process 
in which particle "\" is correlated with particle "2", whereas the second one, represented 
by diagrams diagrams 15c (16c), describes the process in which dynamical correlations are 
acting between particles "1" and '"3". In Figs. El and ^] we show the effect of the hole and 
spectator contributions on the charge density and momentum distributions, respectively. It 
can be seen that the effects on the two quantities are very different: as far as the density 
is concerned, Ap H and Ap s are almost of the same value, and of opposite sign, with a 
small net effect; as for the momentum distribution, the spectator contribution only affects 
the SM distribution by an almost constant factor of 0.8, whereas the hole contribution 
create the large amount of high momentum components. The spectator contribution leads 
to a renormalization of the mean field orbitals and to a decrease of the occupation number 
for states below the Fermi level, whereas the hole contribution is responsible for the high 
momentum components. This explains and qualitatively justifies the parameterized n(k) of 



Ref. 



4l| in the form n(k) = n (k) +ni(k). It is clear, that the amount of hole and spectator 



correlations also depends upon the amount of the mean field contributions; calculations 
show however that the latter, even if obtained within the most sophisticated mean field 
approaches, cannot never provide, e.g. the amount of high momentum components generated 
by the hole contribution, so that the high momentum part of n(k) is practically due only 
to (hole) correlations. Other quantities which are very sensitive to hole correlations will be 



discussed elsewhere 



42|- 



VIII. SUMMARY AND CONCLUSIONS 



In this paper we have addressed the problem of developing a method which could be 
used to calculate scattering processes at medium and high energies within a realistic and 
parameter-free description of nuclear structure. The ^-expansion seems to satisfy such a 
requirement: as a matter of fact, it can be used within the following strategy: i) the values 
of the parameters pertaining to the correlation functions and the mean field wave functions, 
can be obtained from the calculation of the ground-state energy, radius and density of the 
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nucleus, to a given order of the expansion; ii) using these parameters, any scattering process 
can be evaluated at the same order of the cluster expansion. The method therefore appears 
to be a very effective, transparent and parameter-free one. It should however be pointed 
out that, as any other many body approach, our cluster expansion approach may suffer 
from the well known convergence problem, so that the role played by the disregarded higher 
order terms has to be estimated. This is precisely what has been done in the present paper, 
adopting the following procedure: i) our lowest order results have been compared with the 
ones obtained within more complete approaches, like e.g.the FHNC and VMC methods, and 
ii) a direct calculation of the higher order terms of the momentum distribution n(k) has 
been performed. It turned out that the value of the ground-state energy calculated within 
the first order ^-expansion reasonably agrees with the one obtained within the FHNC/SOC 
approach. The agreement is very good as far as the average kinetic energy is concerned, 
whereas differences occur in some of the potential energy contributions, as it should have 
been expected due to the complex spatial dependence of some of the components of the 
nucleon-nucleon interaction. Nonetheless, using the same correlation functions as in the 
FHNC/SOC calculation, we obtain a reasonable minimum value of the energy, with mean 
field WF very near to the ones of the FHNC/SOC approach. Furthermore, our results for the 
charge density and momentum distributions shows a very good agreement with the results 
obtained within the FHNC/SOC approach and even with the VMC approach, and the direct 
calculation of the higher order terms in the expansion of the momentum distributions shows 
a very good convergence of the 7/-expansion up to very high values of the momentum. 

To sum up, we have shown that, using realistic models of the nucleon-nucleon interaction, 
a proper approach based on cluster expansion techniques can produce reliable approxima- 
tions for those diagonal and non diagonal density matrices which appear in various medium 
and high energy scattering processes off nuclei, so that the role of nuclear effects in these 
processes can be reliably estimated without using free parameters to be fitted to the data. 
The approach has already been extended to the treatment of the final state interaction effects 
in A(e, e'p)X processes at medium energies within the eikonal-Glauber multiple scattering 
;heory, and to the calculation of nuclear and color transparencies effects. Preliminary results 



22| are very encouraging. Calculations of other types of high energy scattering processes 



ejr.the total nucleon-Nucleus cross section) are in progress and will be reported elsewhere 
421. 
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APPENDIX A: MEAN FIELD WAVE FUNCTIONS 

The mean field wave functions have the general form: 

= Rni(r) Yije,(p) (Al) 

with R n i{r) the radial part and Yi m (9, if) the spherical harmonics. We have used harmonic 
oscillator and Saxon- Woods wells to generate the radial part; in the harmonic oscillator case, 
we have 

Rni(r) = e-^x^U^X'J 2 ^ (A2) 

with x = r 2 /a 2 , a is the HO parameter and 

» {-l) k x k n\2 k (2Z + 1)!! 

h (n-k)\k\(2l + 2k + l)\\ ' 1 ' 



2'-«+ 2 (2/ + 2n + 1)!! 
(2/ + l)!! 2 n!vr 1 /2 a 3 

whereas, in the Saxon- Woods case, the radial part is the solution of the radial Schrodinger 
equation with one-body potential of the following form 

WO = - 1 + e -l°- Ro)/ao ■ (A5) 

APPENDIX B: PARAMETERIZATION OF THE CORRELATION FUNCTIONS 

FOR 16 AND m Ca 

The correlation functions f^ n \r) for 16 and 40 Ca, shown in Figs. El and 01 respectively, 
can be conveniently parameterized in the following way: 

f {n \r) = E4 n) ^ f>r ' n = l,...,5 (B2a) 



i=i 

/( 6 )(r) = 46) r 2 e -<V + j2 Afr^e-^^ (B26) 

i=2 
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where the parameters A\ and B\ are given in Table HXl for 16 and in Table IXl for 40 Ca. 

APPENDIX C: THE COEFFICIENTS OF THE ONE BODY NON DIAGONAL 
DENSITY MATRIX RESULTING FROM THE SPIN-ISOSPIN SUMMATION 

The coefficients appearing in Eq. ()41|) for the OBMD are defined as 

C%' q) (r 12 ,r V2 ) = £ (<TiT 1 <T2T2\d& ) d$\<r 1 T 1 <r 2 T 2 ), 

C { £' q \r 12 ,r V2 ) = (^i^ldgdgl^tnn) (CI) 

0"1,0-2,T1,T2 

and can be calculated analytically; their explicit values are summarized in Table IXI1 

The coefficients appearing in the definition of the TBD are more involved and can only 
be written in terms of the spin-isospin states upon which they have to be calculated: 

A^ 1>6) = {16, 0, 0, 0, 0, 0} ; AST 1 ' 65 = {4, 12, 12, 36, 0, 0} ; (C2) 

B t$P = E fan, ^r 2 , ^ 3 t- 3 | Off OS? off \a in , a 2 r 2 , a 3 r 3 )) P (C3) 

cr,r 

C(n)!i = E fan, <r 2 r 2 , a 3 r 3 | 6g 6«? 6g K^n, a 2 r 2 , ^3)^ (C4) 

cr,r 

= E fan °"2r 2 , a- 3 r 3 , cr 4 r 4 | O^ O34 Kern, cr 2 r 2 , cr 3 T 3 , 0-4x4) )p (C5) 

cr,r 

where only linked permutations are considered; for example, in the four-body term, the 
identical permutation |ai/?273&i) is not linked, because the only present links are between 
particles 12 and 34, but the two clusters are not connected; there are four unlinked permu- 
tations in this term. 
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TABLE I: The results of the benchmark calculation of the ground-state energy of 16 using the 
V8' interaction 



n 



361 ] , the correlation functions shown in Fig. ^ t_9 ! ] , and harmonic oscillator mean 
field wave functions with parameter a = 2 fm (cf. Appendix Al). The results of the ^-expansion 
obtained in this paper are compared with the FHNC/SOC results of Ref. [9(. (V) is the average 
potential energy, (T) the average kinetic energy, E = (V) + (T) the total energy, and E/A the 
total energy per particle. The kinetic energy of the Center-of-Mass motion has been subtracted 
from the expectation value of the kinetic energy operator. All quantities in MeV . 



approach 


(V) 


(T) 


E 


E/A 


ry-expansion, this paper 


-390.37 


323.50 


-65.90 


-4.12 


FHNC/SOC, Ref. [9] 


-390.30 


325.18 


-65.12 


-4.07 



TABLE II: The contributions (V^) of the first six channels of the V8' interaction to the average 
potential energy (V) shown in Tabled All quantities in MeV. 



approach 


(Vc) 


(Va) 


(Vr) 


(Var) 


(Vt) 


(Vtr) 


v = EiV 


7/-expansion, this paper 


0.6 


-35.4 


-10.1 


-172.8 


-0.03 


-172.7 


-390.37 


FHNC/SOC, Ref. [9J 


0.7 


-40.1 


-10.6 


-180.0 


0.07 


-160.3 


-390.30 
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TABLE III: The results of the benchmark calculation of the ground-state energy of Ca using 
the six-component truncated Urbana UlA interaction, the correlation functions shown in Fig. |2] 
Q], and harmonic oscillator mean field wave functions with HO parameter a = 1.654 /m(see Eq. 
(|A2jO . The results of the //—expansion obtained in this paper are compared with the FHNC/SOC 
results of Ref. |8j. Notations are the same as in Table I. All quantities in MeV. 



approach 


(V) 


(T) 


E 


E/A 


r\ expansion, this paper 


-1655.15 


1425.90 


-229.25 


-5.73 


FHNC/SOC, Ref. [8] 


-1891.60 


1587.60 


-314.80 


-7.87 



TABLE IV: The contributions of the first six channels of the t/14 interaction to the average 
potential energy (V) shown in Table. II 111 Notations are the same as in Table I. All quantities in 
MeV. 



approach 


(Vc) 


(Va) 


(Vr) 


(Var) 


(Vt) 


(Vtr) 


V = Y,iVi 


77-expansion, this paper 


-14.57 


83.20 


91.93 


-1353.45 


11.61 


-473.87 


-1655.15 


FHNC/SOC, Ref. [8] 


-8.40 


92.00 


108.40 


-1549.20 


11.60 


-565.60 


-1891.20 
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TABLE V: The results of the calculation of the ground-state energy and radius of O using the 
full V8' interaction, the correlation functions shown in Fig. |3] [7( and HO and SW mean field wave 
functions. The value of the HO parameter is a = 2.0 fm and the parameters of the SW well are as 
follows: V = 42.0 fm, R = 3.6 fm and a Q = 0.55 fm. The results of the 77-expansion obtained 
in this paper are compared with the FHNC/SOC results of Ref. L 7]. Notations are the same as 
Table I. < r 2 > is the rms radius. Energies in MeV, radii in fm. 



Mean Field 


Approach 


(V) 


(T) 


E 


E/A 


< r 2 > 1/2 


HO 


77-expansion, this paper 


-420.39 


350.39 


-67.54 


-4.40 


2.99 


HO 


FHNC/SOC, Ref. [7] 


-439.84 


353.44 


-86.40 


-5.40 


3.03 



TABLE VI: The same as in Table |V| but for Woods-S axon mean field wave functions. 



Mean Field 


Approach 


(V) 


(T) 


E 


E/A 


< r 2 > 1/2 


SW 


77-expansion, this paper 


-500.59 


444.10 


-56.50 


-3.50 


2.64 


SW 


FHNC/SOC [2 


-519.68 


428.16 


-91.52 


-5.72 


2.83 
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TABLE VII: The same as in Table IVl for 40 (7o; the value of the HO parameter is a = 2.1 fm. 



Mean Field 


Approach 


(V) 


(T) 


E 


E/A 


2 1/2 


HO 


^-expansion, this paper 


-1320.22 


1048.22 


-272.00 


-6.80 


3.72 


HO 


FHNC/SOC [7J 


-1521.20 


1193.60 


-327.60 


-8.19 


3.65 



TABLE VIII: The same as in Table EH but for the Woods-S axon well with parameters V a = 50.0 
fm, R = 5.3 fm, a Q = 0.53 fm. 



Mean Field 


Approach 


(V) 


(T) 


E 


E/A 


2 1/2 


SW 


^-expansion, this paper 


-1293.96 


1018.19 


-275.77 


-7.00 


3.75 


SW 


FHNC/SOC [7] 


-1547.20 


1215.20 


-332.00 


-8.3 


3.66 
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TABLE IX: The values of the parameters appearing in the parametrization given by Eqs. IB2alB2fcl 
of the correlation functions of 16 shown in Fig. 3. 



11 


An) 


An) 
A 2 


An) 
^3 


An) 


A (n) 


An) 
^6 


B (n) 


B (n) 
n 2 


B (n) 


o(n) 

-D4 


B (n) 


1 


1.0005 


0.37314 


-1.1781 


0. 


0. 


0. 


1.0 


2.0 


0. 


0. 


0. 


2 


0. 


-0.1372 


0.1916 


-0.0226 


-0.0141 


0. 


5.0 


3.5 


1.0 


0.13 


0. 


3 


0. 


-0.0795 


0.1271 


-0.0121 


-0.0330 


0. 


5.0 


3.5 


1.5 


0.14 


0. 


4 


0. 


-0.3817 


0.4863 


-0.0535 


-0.0424 


0. 


4.5 


3.7 


1.6 


0.15 


0. 


5 


0. 


0.0114 


0.0527 


-0.0702 


0.0064 


0. 


0.8 


1.5 


1.5 


1.5 


0. 


6 


0. 


-0.1776 


-0.0054 


-0.0237 


-0.00006 


0. 


1.7 


1.0 


1.3 


0.01 


0. 



TABLE X: The same as in Table llXl for the correlation functions of 40 Ca shown in Fig. 4. 
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4 n) 


An) 
A 3 


An) 


An) 
^5 


An) 


B (n) 


B (n) 

D 2 




B (n) 


B (n) 


1 


1.00039 


0.7576 


-1.6015 


0. 


0. 


0. 


3.9 


2.9 


0. 


0. 


0. 


2 


0. 


-0.0573 


0.0965 


-0.0156 


-0.0145 


0. 


7.0 


3.5 


0.8 


0.2 


0. 


3 


0. 


-0.0207 


0.0474 


-0.0019 


-0.0279 


0. 


8.5 


3.5 


2.0 


0.215 


0. 


4 


0. 


-0.0290 


0.0906 


-0.0237 


-0.0356 


0. 


9.4 


3.0 


1.0 


0.22 


0. 


5 


0. 


0.0165 


0.0061 
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TABLE XI: The value of C%' q) (r 12 , r V2 ) and Cf q > (r 12 ,r V2 ) defined in Appendix C. The order 
of the operator p, q = 1, 2, . . . , 6 is the same as Table UTI Here (Si 2 S\/ 2 ) is defined as (Si 2 Si/ 2 ) = 
12(3(f 12 ■ fi'2) 2 — 1), with the definition of r = r/r. 
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? IG. 1: The correlation functions for 16 (9 correspondine to the truncated Argonne AV8' interaction 

i n 

36l | used in the benchmark calculation (After Ref. |9|]). 
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^IG. 2: The correlation functions for 40 Ca corresponding, to the truncated Urbana U14 interaction 
3tI ] used in the benchmark calculation 36] (After Ref. ]{]]). 
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FIG. 3: The correlation functions for 16 corresponding to the ^4^8' interaction (After Ref. 
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FIG. 4: The correlation functions for 40 Ca from .7], corresponding to the AV8' interaction 
(After Ref. Q). 
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FIG. 5: The charge density of ie O. Thick full: experimental density. Thin full: results of the 
7/-expansion with harmonic oscillator (HO, top) and Saxon- Woods (SW, bottom) wave functions 
and correlations functions shown in Fig. |31 The wave function parameters correspond to the 
minimization of the ground-state energy. Dots: mean field density obtained by setting /W = l ; 
j-{n^i) _ q rp^g chaj-ge density is obtained by folding the matter density with the charge density 
of the proton and correcting for the center-of-mass motion effects. The value of the rms radius is 
(r 2 )V 2 _ 3 07 f m , with HO wave functions, and (r 2 ) 1 ^ 2 = 2.85 fm, with SW wave functions. The 
value of the HO parameter is a = 2.00 fm and the parameters of the SW well are V Q = 42.0 fm, 
R = 3.6 fm and a Q = 0.55 fm. The density normalization is J drp(r) = Z. 
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FIG. 6: The same as in Fig. but for 40 Ca, and correlation functions from Fig. |1J The value of 
the rms radius is (r 2 ) 1 ^ 2 = 3.72 fm, with HO wave functions, and (r 2 ) 1 ^ 2 = 3.75 fm, with SW 
wave functions; the value of the HO parameter is a = 2.10 fm, and the parameters of the SW well 
are V Q = 50.0 fm, R a = 5.3 fm and a Q = 0.53 fm. 
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FIG. 7: The ground-state energy of O versus the harmonic oscillator parameter a calculated with 
the ^-expansion and the V8' interaction, using the correlation functions shown in Fig.3. 
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FIG. 8: The same as in Fig, 03 but with mean field wave functions chosen so as to better reproduce 
the experimental density. The value of the rms radius is (r 2 ) 1 ^ 2 = 2.73 fm, with HO wave functions, 

1 /2 

and (r 2 ) = 2.71 fm, with SW wave functions. The value of the HO parameter is a = 1.81 fm, 
and the parameters of the SW well are V a = 53.0 fm, R a = 3.45 fm and a Q = 0.7 fm. 
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FIG. 9: The same as in Fig. El but for 40 Ca. The value of the rms radius is (r 2 ) = 3.56 fm, 
with HO wave functions, and (r 2 ) 1//2 = 3.34 fm, with SW wave functions. The value of the HO 
parameter is o = 2.00 fm and the parameters of the SW well are V = 50.0 fm, R = 5. fm and 
a a = 0.515 fm. 
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FIG. 10: The two-body densities (Eq. (j42jl ) of 16 (9 and 40 Co corresponding to the correlation 
functions of Fig. |3] ( ie O) and Fig. 21( 40 Ca); the quantities in the figure are integrated over the 
center of mass coordinate (see Eq. (|51|0 , The splitting of the two-body density in n different 
quantities is explained in Eq. (|42|) and the following text. 
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FIG. 11: The momentum distributions of 16 corresponding to harmonic oscillator (top) and 
Saxon- Woods (bottom) wave functions giving the best density shown in Fig. |HJ The full thin 
line includes only the central correlation function, whereas the thick full line includes all of them. 
Our results are compared with the results of Ref. 3] (stars), obtained with the same correlation 
functions. The results of Ref. [5] obtained within the VMC approach using the AV14 interaction 
are also shown by full squares. The value of the kinetic energy obtained by integrating n(k) are: 
(T) = 297.87 (central, HO), (T) = 476.55 MeV (full, HO); (T) = 306.99 (central, SW) and (T) = 
494.48 MeV (full, WS). The normalization of n(k) is / dkn(k) = 1. 



38 




FIG. 12: The same as in Fig. 1111 but for 40 Ca and correlation functions from Fig. 0]and mean field 
wave functions giving the best charge density of Fig. |5J The value of the kinetic energy obtained 
by integrating n(k) are (T) = 782.87 (central, HO), (T) = 1178.45 MeV (full, HO); (T) = 836.24 
(central, SW) and (T) = 1245.21 MeV (full, SW). 
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FIG. 13: The effect of the various correlation functions on the momentum distribution of O. fi 
approximation: only central correlation; f% approximation: = = = 0. approxi- 

mation: full correlation set, n = 1,...,6. Calculations were performed with correlation functions 
taken from Fig. 3 and HO wave functions. 




FIG. 14: The convergence of the momentum distributions of O calculated by considering the 1st 
and 2nd orders of the ^-expansion, Saxon- Woods mean field wave functions, the f% approximation 
and the correlation functions of Fig. |31 Our results are compared with the results of Ref. [j| 
obtained within the VMC approach and the AV14 interaction. The values of the kinetic energies 
obtained by integrating the momentum distributions are: (T) = 521.87 MeV (mean field), (T) = 
980.10 MeV (1st order rj), (T) = 932.64 MeV (2nd order rj). The normalization of n(k) is 
Jdkn(k) = 1. 
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FIG. 15: Diagrammatic representation of the one body density p(ri) in the lowest order of the 
77-expansion (Eq. (|5t)j)). Dots denote the corresponding spatial coordinates and integrations occur 
over full dots. An oriented full line represents the shell model uncorrelated non diagonal density 
matrix p<?\i,j), an oriented closed line the diagonal density matrix p (i,i), and a dotted line the 
correlation H(rij,rji) between particles i and j. a): shell model, uncorrelated density p Q ; b): hole 
contribution Ap H c): spectator contributions Ap s . The direct and exchange contributions are 
shown on the left and right sides of the Figure, respectively. 




FIG. 16: The same as in Fig. EDf° r the one body mixed density matrix p^(r,r'). 
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FIG. 17: The charge density of 16 (9. Thin full: shell model density (Eq. I| 33|) ). Thick full: 
correlated density (Eq. (|50j)) calculated with the correlation functions of Fig. El and HO mean 
field Wave Functions with parameter a = 1.8 fm. Dashed: hole contribution Ap H (r); dotted: 
spectator contribution Ap (r). 




FIG. 18: The same as in Fig. 1171 but for the momentum distribution. 
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